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ABSTRACT 

Context. Using Herschel data from the Open Time Key Project the Herschel Virgo Cluster Survey (HeViCS), we investigated the 
relationship between the metallicity gradients expressed by metal abundances in the gas phase as traced by the chemical composition 
of HII regions, and in the solid phase, as traced by the dust-to-gas mass ratio. 

Aims. We derived the radial gradient of the dust-to-gas mass ratio for all galaxies observed by HeViCS whose metallicity gradients 
are available in the literature. They are all late type Sbc galaxies, namely NGC4254, NGC4303, NGC4321, and NGC4501. 
Methods. We fitted PACS and SPIRE observations with a single-temperature modified blackbody, inferred the dust mass, and calcu- 
lated two dimensional maps of the dust-to-gas mass ratio, with the total mass of gas from available HI and CO maps. HI moment- 1 
maps were used to derive the geometric parameters of the galaxies and extract the radial profiles. We examined different dependen- 
cies on metallicity of the CO-to-H2 conversion factor (Xco), used to transform the '^CO observations into the amount of molecular 
hydrogen. 

Results. We found that in these galaxies the dust-to-gas mass ratio radial profile is extremely sensitive to choice of the Aco value, since 
the molecular gas is the dominant component in the inner parts. We found that for three galaxies of our sample, namely NGC4254, 
NGC4321, and NGC4501, the slopes of the oxygen and of the dust-to-gas radial gradients agree up to ~O.6-O.7R25 using Aco values 
in the range 1/3-1/2 Galactic Aco ■ For NGC4303 a lower value of Acq ~0. 1 x 10^" is necessary. 

Conclusions. We suggest that such low Acq values might be due to a metallicity dependence of Acq (from close to linear for 
NGC4254, NGC4321, and NGC4501 to superlineai' for NGC4303), especially in the radial regions Rq <0.6-0.7R25 where the molec- 
ular gas dominates. On the other hand, the outer regions, where the atomic gas component is dominant, are less affected by the choice 
of Aco , and thus we cannot put constraints on its value. 

Key words. Galaxies: individual: NGC4254, NGC4303, NGC4321, NGC4501 — Galaxies: spiral — ISM:dust, extinction — 
Galaxies: ISM - Galaxies: abundances — Submillimeter: galaxies 



1. Introduction 



Virgo is one of the best studied galaxy clusters, being the rich- 
est cl uster nearest to our own Galaxy (-17 Mpc. iGavazzi et al] 
I1999 I1. It is a relatively populous system, consisting of more 
than 1000 confirmed members (Binggeli et al. 1985). Galaxies 
in clusters such as Virgo differ significantly from their field 
counterparts since interactions with the hostile environment re- 



move gas, quenching the star formation (c.f.. lBoselli & Gavazzil 
l2006l) . 

A galaxy's metallicity is closely related to the star formation 
(SF) history by which the interstellar medium (ISM) is enriched 
with the end-products of stellar evolution, and to the infall pro- 
cess that dilutes the ISM and triggers new SF. As a consequence 
of their star formation histories, gas stripping and infall events, 
modified by the cluster environment, galaxies in clusters are ex- 
pected also to diff'er in metal content relatively to isolated galax- 
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ies. A fundamental tool with which tracing the chemical evo- 
lution of a galaxy is the study of its radial metallicity gradient. 
The metallicity gradient tracks indeed the star formation history 
of galaxies, integrated over time, together with infall and/or out- 
flow events. 

The first pionee ring work in Virgo was done by 
ISkillman et al. I (Il996l) . who analyzed nine spiral galaxies with 
the aim of seeking correlations among their gas content, lo- 
cations in the cluster, metallicities an d rad ial gradients, and 
comparing them with field spirals. ISkillman et al. I (11996) found 
weak evidence of shallower gradients in cluster galaxies defi- 
cient in HI than gradients in galaxies with a normal HI content. 
The situation is however very complex because galaxy interac- 
tions affect the star formation history and the gas content across 
the disk, producing metallicity gradients which differ from those 
measured in isolated galaxies of the same morphological type. 
iRupkeetal. I (120101) have shown that gradients in strongly inter- 
acting galaxies are flatter than in similar isolated galaxies; on the 
other hand, the metallicity gradient of M8 1 , the largest member 
of a small group of galaxies, is steeper than i n an isolated coun- 
terpar t due to gas removal in the outskirts (IStanghellini et al.l 
l2010h . 

A direct correlation between gas metallicity and the dust-to- 
gas mass ratio is naturally expected since approximatively half 
of the metals in the ISM reside in dust grains; thus the dust-to- 
gas ratio should scale with metal abundance. Such a trend has 
been obtained theoretically by models computing consistently 
the evolution of metals and dust, desp i te the large u ncertain- 
ties in the yields of both (Dwek' 1998; 'lnouel l2003h . The re- 
lation of the global dust-to-gas mass rati o with metalhcity wa s 



investigated by, e.g., James et al. J 2002h; Draine et alTI (l2007h ; 



iHirashita et aD (l2008h ; iLisenfeld & Ferrara I (Il998l) . The radial 
va riation of the dus t-to-gas mass ratio was first investigated 
bv 'Issa et al.1 (11990) in our Galaxy and in other nearby galax- 
ies (LMC, SMC, M31, M33, and M51) who found evidence 
for a correlation, with dust-to-gas mass ratio and metallic- 
ity decreasing at roughly the sa me rate with inc reasing galac- 
toce ntric radius. More r ecently, iBoissier et al. I (j2004, 2005) 
and iThilker et al. I (l2007h found a clear relationship between 
metallicity and extinction, thus dust, in several nearby galax- 
ies, suggesting that the va riation in extinction i s asso ciated 
with the metallicity gradient. lMufioz-Mateos et aL 1 (l2009l) found 
a good correlation between dust-to-gas and metallicity gradi- 
ents in the Spitzer In frared Nearby Galaxie s Surve y (SINGS; 
iKennicutt et al. Il2003h . Finally, Bendo et al. I ( l2010al) compared 
the dust-to-gas ratio and metallicity gradients in NGC2403, find- 
ing a similar decreasing behavior with radius. 

We have recently obtained observations of the Virgo galax- 
ies with the Herschel Space Observatory (Pilbratt et al. 2010), 
within the Ope n Time Key Projec t HeViCS (Herschel Virgo 
Cluster Survey) (iDavies et al. Il2010b . HeViCS maps a wide area 
over the Virgo Cluster at wavelengths from 100 to 500 fim. This 
spectral range covers the peak of the thermal emission from cold 
dust (T<30 K) which enables the detection of the bulk of the dust 
emission in galaxies. Also, Herschel gives an unprecedented res- 
olution at these wavelengths (ranging from about 10 to 36", 
equivalent to 1-3 kpc for Virgo galaxies). Atomic and molec- 
ular gas maps are available in the literature at a comparable res- 
olution (Chung et al. 2009a; Kuno e t al. .2007.) . thus providing 
resolved maps and the possibility of deriving radial gradients of 
dust-to-gas mass ratios. 



In the present paper we assessed the validity of using 
HeViCS observations to obtain metallicity gradients from radial 
profiles of the dust-to-gas mass ratio. We investigated the hy- 
pothesis that the local dust-to-gas mass ratio is proportional to 
metallicity, starting from the four spiral galaxy in the Virgo clus- 
ter (NGC4254, NGC4303, NGC4321, and NGC4501) whose 
metallicity gradients are available in the literature. We studied 
the relation between radial profiles of metallicity and dust-to-gas 
mass ratios, and how this can be used to constrain the CO-to-H2 
conversion factor (Xco ) and its dependence on metallicity. 

The paper is structured as follows: in Sect. 2 we describe 
the new HeViCS observations, and the observations of atomic 
and molecular gas. We also derive the oxygen abundance and 
its radial gradient in each galaxy. Sect. 3 discusses the fits to 
the observed Herschel dust spectral energy distributions (SEDs), 
how the radial profiles are derived, and how we calculate dust 
and gas masses. In Sect. 4 we develop the method by which we 
constrain Xco , and Sect. 5 gives our conclusions. 



2. Data 

2. 1 . The sample 

In the HeViCS field there are four galaxies for which the 
oxyg en gradient has been well determined in th e literature 
(e.g., ISkilknaneta"nil996t iMoustakas etaTI 120101) . They are 
NGC4254 (M99), NGC4303 (M61), NGC432I (Ml 00), and 
NGC4501 (M88), all Sbc late-type galaxies. 

NGC4254 is a bright spiral galaxy located at the periph- 
ery of the Virgo cluster, at a projected distance of ~1 Mpc 
from the cluster center. Optical images show that this galaxy 
has one-armed structure, also seen in the HI gas distribution 
dPhookun et al. II1993I; IChung et al. "2009a). Such an asymmet- 
ric spiral pattern is often observed in tidally galaxies , but there is 
no apparently massive companion near NGC4254 (ISofue et al. I 
2003.) . 

NGC4321 is located at a distance of -1.1 Mpc form M87, and 
has an HI disk that is slightly larger than the optical disk. 
NGC4501 is the closest galaxy of our sample to M87, being 
located at a distance of ~ 0.5 Mpc. It i s weakl y Hl-deficient, 
following the definition of 'Chung et al. I (l2009a !). Comparisons 
with simulations suggest that NGC45 01 is in an early stage of 
ram pressure stripping (iVolkner et al.1l200"8l) . entering the high- 
density region of the cluster for the first time. 
NGC4303 is the most isolated galaxy in our sample. It is a barred 
spiral galaxy with face-on geometry located in the outskirts of 
the Virgo cluster We ass ume for all the galaxies the distance of 
-17 Mpc (iGavazzi et al.lfr999b . 

2.2. Herschel observations 

The HeViCS program consists of Herschel observations of an 
area of about 60 sq. deg. over the denser parts of the Virgo 
Cluster. The total area is made of 4 overlapping fields, which are 
observed in parallel scanning mode (fast scan rate: 60"/s) with 
both the PACS and SPIRE instruments, yielding data simulta- 
neously in 5 spectral bands, at 100 and 160 fim (from PACS) 
and at 250, 350, and 500 /jm (from SPIRE). At the comple- 
tion of the program, each field will be covered with 8 scans 
done in two perpendicular scan directions. The full width half 
maximum (FWHM) of the beams is 6'.'98 x 12'.'7 and 11'.'64 
X 15'.'65 in the two PACS bands (PACS Observers' Manual, 
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2010), and 18'.'2, 24'.'9, 36'.'3 in the three SPIRE bands (SPIRE 
Observers' Manual, 2010). At the time of writing, each field has 
been observed with at least two scans. In this paper, we use this 
dataset, whose data reduction and an alysis is described in details 
in PAPER VIII. lDavies etaTl JIoTlh . 

Some papers presenting Herschel observation s for these 
galaxies have been a lready published. Eales et al~l (1201 Oh and 
ISauvage et aO (l2010l) analyzed the SPIRE maps of NGC4254 
(M99) and NGC4 32I (MlOO) o bserved within the Herschel 
Reference Survey (iBoselli et al. 112010 ) to map the ISM using 
dust e mission. Adding thes e data to archival Spitzer, HI, and CO 
maps. 'P ohlen et al~1 ( 20101) investigated the spatial distribution 
of gas and dust in these same galaxies. They also present as 
a preliminary result, the ratio of the total gas mass (HI + H2) 
to 500 fim flux, an approximation of the dust mass for the two 
galaxies. They found a decreasing dust-to-gas mass ratio with 
radius, consistent with results by, e.g., iBendo et al. I (I2010al) in 
NGC2403. With the present- time availability of the PACS data 
the dust SED fitting con be better defined allowing to measure 
the exact shape of the radial dust-to-gas mass gradient. Finally, 
[Smith et al. (2010) presented a resolved dust analysis of three 
of the largest (in angular size) spiral galaxies in HeViCS, among 
them NGC4501. 



2.3. The calibration of metallicity and tlie abundance 
gradients 

The metallicity measurements of these galaxies are available in 
the literature from optical spectroscopy of their HII regions. The 
most direct method to derive the oxygen abundance is to measure 
the electron temperature (T^) of the ionized gas using the inten- 
sity (relative to a hydrogen recombination line) of one or more 
temperature-sensitive auroral lines such as [O III] /14363A, [N 
II U 5755A, [S III1 /16312A , and [O II] ^7325 A, as summarized 
bv lMoustakas et ain(l2010t) . However, measurements of the elec- 
tron temperature were not available fr om th e spectroscopic 

(198j 



observations in the original paper s of i McCal l et al. , _ . 
[Shields etaD ( 1199 Ih . iHenrv et al. I (1 19941) . and ISkillman et al. 
(il996) . The temperature diagnostic lines are indeed intrinsically 
faint in metal-rich HII regions. Therefore, oxygen abundance 
has been derived using the strong-line abundance calibrations 
which relate the metallicity to one or more line ratios involv- 
ing the strongest recombination and forbidden lines. In particu- 
lar, the oxygen excitation index R23=([OII]h-[OIII])/H^ is one of 

the most often adopted calibrators to estimate the nebular abun- 

I 1 ' ' 

dances. As discussed bv Moust akas et al. ' (2010), the principal 

advantage of R23 as an oxygen abundance diagnostic is that it is 
directly proportional to both principal ionization states of oxy- 
gen, whereas one of the major disadvantages is that the rela- 
tion between R23 and metallicity is degenerate for low and high 
metallicity. 

However, the calibrations of the metallicity by means 
of the strong-line ratios are not unique. They can be di- 
vided in three main categories: those calibrated with photoion- 
ization models (e . g., McGaugh 1991; Zaritsky et al. 199^ 
iKewlev & Dopital llOQ'A iKobulnickv & KewlevI 120041 : those 
calibrated directly wit h the electron t e mperature, called em 



scale. The oxygen abundances derived using the theoretical cal- 
ibration are up to a factor of ~4 higher than those based on 
the empirical calibration. However, despite the significant zero- 
point offset in the abundance scales, to first order the slope of the 
abundance gradients agree when c alculated with different cali- 
brators ( see lMoustakas et al. ]l2010l) . 

We have tested the effect of several metallicity calibrators in 
NGC4254, the galaxy with the best sampled m etallicity gradient. 
Starting with the abundance estimates from iMoustakas et al. I 
(2010), based on literature spectroscopy calibrated with the 
Kobulnicky & Kewley (2004) form ula (hereafte r KK04), we 
used the relationships provided by Kewlev & Ellison I (l2008h 
to convert to other abu ndance scales. The relationships of 
iKewlev & Ellison 1 (l2008l) are obtained for limited oxygen abun- 
dance ranges, corresponding to the ranges where they could per- 
form a polynomial fit to transform one abundance scale to an- 
other Because of this, from the abundances of KK04 we were 
unable to r ecover the oxygen d etermination of iPettini & Pagel I 
(l2004l) and iPilvugin & Thuanl (12005). For these cases, we re- 
compute d the oxy gen abundanc e from the original spe ctra 
(McCall etal. 111985 : Shields et a Ol99ll:lHe"nrv et al. 111994 . 

In Fig. [H we show the result of our test: while the slope of 
the gradient is almost invariant with different calibrations, the 
zero-point depends on the choice of the cjilibration. In partic- 
ular, the empirical calibra tions of iPettini & Pagel I (12004 ) and 
iPilvugin & Thuan ] (l200l show values of the oxygen abun- 
dances roughly 0.5-0.8 dex lower than the other determinations. 
These two latter empirical calibrations were obtained for HII 
regions with available electron temperature in a relatively low 
metallicity regime. They could not be valid for metal rich en- 
vironments, as the galaxies of our sample. In Sect. 14.1 1 we will 
discuss how the dust-to-gas ratio might help in setting a lower 
limit to the metallicity and to discriminate among different cali- 
brations. 



NGC4254 




Fig. 1. The oxygen abundance gradient of NGC4254 obtained 
with several metallicity calibrations: Kobulnicky & KewlevI 



pirical mefliods (e g.. IPilvugin 1 120011: IPettini & Pagel] 12004 (l2004 ( magenta circles ) [Pilwgin & Thuanl (l2005l) (veUow 
trr; ■■ ^ ^^""'^ ■ squares), IZaritskv et al. I ( 1994*) (cya n trian gles), iMcGaughl 



IPilvugin & Th uan' '2005^ 



ods (e.g.. 



and those combining both meth 
2002). As discussed widely in 



j Denicolo et a l 

IKewlev & Ellison" T (t2008h . the abundances derived with differ 
ent methods do not have a common absolute oxygen abundance 



( 199 11) (blue stars). lKewl"ev& Dopita] (|2002|) (empty squares) 
Tremont i et al. ] (12004 (green asterisks), and IPettini & Pagel I 



(2004) (violet diamond). 
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2.4. HI and CO maps 

The radial profile of the gas, including both atomic and molec- 
ular components, is necessary to derive the dust-to-gas mass ra- 
tio gradient. For NGC4254, NGC4321, and NGC4501, we use 
the moment-0 HI maps obtained with VLA Imagi ng survey of 
Virgo g alaxies in Atomic gas (VIVA) survey by Chung et al. I 
(l2009ah . VIVA observations reach a column density sensitiv- 



ity of 3-5x10 cm . The comparison of their total HI fluxes 
with values in the literature from single dish observations gives 
a good agreement especially for the large galaxies, indicating 
n o loss of flux in the interferometric observations (see Fig. 5 
in lChung et alH Tooga). The beam sizes are: 3T86x 28'.'07 for 
NGC4254, 15'.'90x 14'.'66 forNGC4321, and 16'.'83x 16'.'41 for 
NGC4501. The HI ra dial profile of NGC4303 is available from 
IWarmels I (1 19881) and lCavatte et al. I (fT990l) . We adopt flie com- 
bined radial profile of the two, shown in Fig. 3 of ISkillman et al. I 
(119961) . 

Maps of molecular gas were available thanks to the '^CO (J 
= 1-0) mapping survey of 40 nearby s piral galaxies, perf ormed 
with the Nobeyama 45 m telescope bv lKuno et al. I (l2007h . 



3. Analysis 

3.1. The mass and temperature of dust 

Maps of the dust temperature a nd mass surface density were 
obtained as in ISmith et al. I IMO ). The images of the galaxies 
from the five PACS (100 and 160 yum) and SPIRE (250, 350, 
and 500 fim) bands were all convolved and re-gridded to the 
lower resolution (FWHM=36'.'9, 3 kpc) and pixel size (14'.'0, 
1.1 kpc) of the 500 fim observations. We used only pixels with 
S/N> 10 at 50 ^m (the rms of ou r maps in this band is about 
0.3 MJy/sr, see iDavies et al. 1 120111) . The selection of these high 
surface-brightness pixels was necessary to limit the uncertainties 
due to background subtraction and avoid the artefacts caused by 
the high-pass filtering in the PACS data reduction. Despite this 
limit, we were able to study the dust and gas properties up to 
at least 0.7 in all galaxies. For each galaxy we thus con- 
sidered approximately 150-200 pixels covering about 1/2 of the 
optically defined area. 

We estimated the error on the surface brightness on a pixel 
by pixel basis by comparing galaxy images from the two scan 
data used in this paper with those relative to other two scans 
recently taken by Herschel which cover only part of the He ViCS 
field. We f ound errors very sim ilar to those estimated on the total 
fluxes (see iDavies et al. I I20T1I) . Including a calibration error of 
15% for PACS (PACS ICC, private communication) and 7% for 
SPIRE (SPIRE Observer' Manual, 2010), the total error is 30%, 
20%, 10%, 10%, and 15% of the flux at 100, 160, 250, 350, and 
500 jum, respectively. 

The SED for each pixel was fitted with a single modified 
blackbody, using a power law dust emissivity k^i - KoiAo/A'f, 
with spectral index /3 - 2 and emissivity kq - 0.192 kg"' at 
A() = 350 jum. These values reproduce the averag e emissivity o f 
models of the Milky Way dust in the FIR-submm ( iDraine I2003h . 
The fit was obtained with a standard minimization technique. 
In the pipeline calibration, the flux density observed by the var- 



' R25 is the radius of the galaxy measured to a B surface brightness 
of 25 mag arcsec"^, and is an indication of the size of the galaxy; R25 



ious instruments, i.e. weighted over each filter passband, is con- 
verted into a monochromatic flux density assuming Fy cc y^' 
Before fitting a modified blackbody, a color correction should 
be applied to the data, to account for the real spectral slope of 
the source. Alternatively, the conversion implemented into the 
pipeline calibration can be removed from the data (in SPIRE 
parlance this is equivalent to dividing the pipeline flux densi- 
ties by the K4 factor; SPIRE Obsei-ver' Manual, 2010); the pass- 
band weighted flux thus obtained should then be compared with 
the mean of the model flux density over the spectral response 
function for each of the bands. We adopted this second tech- 
nique, using the appropriate response functions for the PACS 
and SPIRE bands (for SPIRE, we used the response functions 
for extended emission). However, color corrections are small for 
the adopted emissivity and the temperature range derived here 
(see .Davies et al. , 201 1.) . 




100 



500 



X [ iim ] 



were obtained from NED; 
our Galaxy. 



"O.7R25 is equivalent to the solar radius in 



Fig. 2. Typical SEDs on three different positions (pixels) on 
NGC4254. Blue dots are the measured (color-corrected) fluxes 
and the red (solid) + curves are the modified blackbody fits. The 
three SED correspond to minimum 21K), mean 23.5K), 
and maximum (» 26K) temperature in the galaxy (see Fig. 0. 
The green (dashed) curves show the fit of the.Draine & Li (20()7|) 
model to the data, for IRSFs of intensities 2, 4 and 8 times the 
Local ISRF (see their paper for details). 



In Fig.|2]we show typical SEDs and graybody fits for three 
positions of different dust temperature in NGC4254. As shown 
by the figure, a single temperature modified blackbody with 
/5 - 2 is sufficient to obtain reasonable fits of the SED over th e 
wavelength range considered here (see also IDavies et al. |[201 ih . 
When data at shorter wavelength than PACS 100 fim is avail- 
able, one might want to consider a two-component model, to 
include emission from warmer dust that might significantly con- 
tribute at least to the 100 u m flux. This was done, for exam- 
ple, by Bendo et al. I (l2010bh . who used Herschel data, and by 
ISmiflietaL I (120101) who used 70 yum data from the Spitzer satel- 
lite. However, they found that the inclusion of a warmer temper- 
ature component, thought necessary to fit the 70 //m data, im- 
proves the fit at /I > 100 fim only slightly, and does not modify 
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significantly the estimate of the temperature and mass surface 
density of cold dust. 

3.2. Sources of uncertainty in mass and temperature of dust 

When only the errors on photometry are considered, the uncer- 
tainty in the determination of the temperature is about 2 K, and 
thus ~2Q% on the dust mass surface density. 

In principle, fitting the dust SED with a single thermal- 
equilibrium temperature component could result in larger uncer- 
tainties in the dust mass estimates: grains of a given size and 
material could be exposed to different intensities of the interstel- 
lar radiation field (ISRF) and thus attain different equilibrium 
temperatures which will contribute differently to the SED; con- 
versely, for the same radiation field the SED could depend on 
the dust distribution, because it results from the emission of a 
mixture of grains of different size and composition, each with its 
own equilibrium temperature. We found neither of these to have 
a strong effects on our mass estimates. In fact, by fitting the SED 
pixel by pixel, we akeady take into account the gradients due to 
the diffuse ISRF, which is more prominent in the radial direction 
than in the ver tical directions (i.e. a long the line of sight, for non 
edge-on disks: iBianchi et al.]l2000l) . 

However, the temperature radial gradients found here (see 
Sect. 13. 3b would not produce very large uncertainties, even when 
the global SED is fitted with a single temperature model. For 
our targets, the differences between the sum of the dust mass 
in each pixel, and the dust mass obtained by fitting the sum of 
the flux densities in each pixel, is smaller than the fit error The 
insensitivity of the global dust SED fitting on the shallow dif- 
fuse ISRF gradient is clearly shown in the analysi s of the th e 
FIR/submm SED of late type galaxies of Draine e t al. I (l2007h : 
though the more complex fitting procedure include s a dust grain 
model and a range of ISRF intensities (see also iDraine & Lil 
l2007h . the global SED at ^ > 100//m is found to be fitted by 
a dust component that accounts for most of the dust mass (a part 
from a few percents), heated by an IRSF of constant intensity. 

We evaluate d the effects of a co mplex dust mixture by fit- 
ting the m odel of Draine & Lil ( 20071) to our pixel-by-pixel SED. 
Following iDraine et al. I (12007 ). we used a single value of the 
ISRF for emission at A > lOOyum. An example can be seen in 
Fig. [2] The SED from the dust grain mixtures and our single 
temperature model fit equally we ll the data. The du st mass ob- 
tained by using the procedure of IDraine et al. 1 (l2007.) is higher 
by about 10%, that is, within the error we quoted. Thus, the dust 
mass derived from a simple averaged emissivity, and a single 
IRSF - or temperature - , is not severely underestimated. 

Larger uncertainties can come also from the assumption for 
th e emissivity. For exa mple, the value of the emissivity derived 
by I James et al. I (l2002l) from SCUBA observations of galaxies is 
equivalent to kq = 0.41 kg"', a factor of two larger than the 
value we adopted. Adopting this emissivity would result in dust 
surface densities a factor of two smaller than what we found in 
this paper. 

Also, dust emissivity has been reported to increase by about 
a fa ctor two in grains as sociated to denser environments (see, 
e.g.. iBianchi et alJl2003h . However, this might be the case for 
extreme environments and not representative of the bulk of the 
diffuse dust mass: in the Milky Way, recent results from the 
Planck satellite show no emissivity variation betweed dust asso- 
ciated with HI and CO emi ssion, nor emissivity variations w ith 
the galactocentric radius (Planck Collaboration et al.ll201 ll) . In 




Fig. 3. Radial variation of dust temperature in NGC4254 (red 
dot-dashed cui-ve), NGC4303 (green dashed line), NGC4321 
(blue dotted curve), and NGC4501 (black solid curve). 



Sect. 14.1! we will discuss how the uncertainty on kq might affect 
our discussion. 

The limited wavelength coverage does not allow us to inves- 
tigate in details the effect of variation of the emissivity spectral 
index. In any case, the modified /? = 2 blackbody provides good 
fits for all the SEDs analised here, as well as for the globa l SEDs 
in a larger sample of HeViCS objects (iDavies et al. 1201 H) . Local 
variations o f the dust emi ssivity index with temperature, as those 
reported bv lParadis et al.l(i2010l) cannot be easily veryfied in our 
dataset. Variation of f5 from 1 .7 to 2.2, as those measured at a ref- 
erence temperature of 20K, would result in an underestimate and 
overestimate of the dust mass of the same order as the quoted 
errors, respectively. Lower /? values, as those found for dust at 
higher temperatures, will not be able to provide good fits to our 
SEDs in the central part of the galaxies. 

Maps of the temperature and dust mass sur f ace de nsity for 
NGC4501 have been presented in lSmith etaTI (1201 0'). For the 
other galaxies, they will be presented in Vlahakis et al. (201 1, in 
preparation). 

3.3. Radial profiles 

The galaxies in our sample are disturbed by tidal interactions, 
thus we consider more reliable the determination of the geo- 
metric parameters of the galactic disks from kinematics rather 
than from photometry (see, for example, the discrepancy be- 
tween parame ters derived with the two methods for the case of 
NGC4321 in Chung et al. "2009a'). The task ROTCUR of the 
package NEMO (Teuben 1995) was used to fit a tilted ring 
model to the 2 1 -cm moment- 1 maps from the VIVA dataset. We 
obtained values for the rotational velocity, the inclination and the 
position angle, as a function of the galactocentric radius. 

The radial variations of / and p. a. are shown in Fig. |4| The 
average values are indicated with solid lines. For NGC4321 and 
NGC4501, / and p. a. are quite constant with radius and thus 
we use the average parameters in the whole radial range. For 
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NGC4254, which is a one-armed spiral, the variations of p. a. and 
; are important due to its asymmetric shape. Since in our compar- 
ison between metallicity and dust-to-gas mass ratio gradients we 
are interested mostly in the inner regions where the metallicity 
data are available, for NGC4254 we used / and p. a. averaged in 
the regions with galactocentric radius Rc < R25 (da shed Une in 
Fig, m top panels) where the variations are small. Pohle n et al. I 
(I2OI O) also found different ellipse parameters for the outer parts 
and for the inner parts of NGC4254. For NGC4303, for which 
VIVA maps are not available, we adopted the parameters by 
IKoda & Sofuel dloo^ . 

The final values are shown in Table [1] in column 1 we show 
the galaxy name, in columns 2 and 3 the mean / and p. a., in col- 
umn 4 the optical diameter Dgp, in arcmin, in cols. 6, 7, and 8 the 
slope of the metallicity gradient, the metallicity at the equivalent 
solar radius O.7R25, and their reference. For NGC4254 we report 
the value we adopt, i.e. the average within the optical radius, 
and in square brackets the average in the whole radial range. 
Between square brackets, for all galaxies, we report also inclina- 
tions and position angles from the literature, which are generally 
in good agreement with our values in the radial range consid- 
ered. 

The atomic and molecular gas maps were convolved to the 
same resolution of the SPIRE 500 //m map and regridded on the 
same pixel scale. The radial bin size in the profiles correspond to 
1 pixel in the maps (14" corresponding to ~1.1 kpc). Adopting 
the geometric parameters in Table [1] we used a dedicated IDL 
procedure to obtain radial profiles (azimuthally averaged over 
elliptical annuli) for HI, CO, and dust mass surface density. The 
profiles are shown in Figs. |5] the top panels contain PACS and 
SPIRE profiles at 100, 160, 250, 350, and 500 yum (from top to 
bottom) for NGC4254, NGC4303, NGC4321, and NGC4501; 
the bottom panels the profiles of dust (green), HI (red), and H2 
(blue) surface density, where the CO flux is converted into H2 
using the standard conversion factor Xco~ (1-8 ± 0.3) x 10^" 
cm"^ (K km s"')"' (see below). For all galaxies H2 is dominant 
in the inner-most regions, while HI becomes dominant over H2 
at R>0.5-0.6 R25. HI radial profiles are quite flat within R25, 
whereas both H2 and dust profiles are decreasing with radius. 

The errors on the radial profiles were estimated by a com- 
bination in quadrature of the uncertainty in the overall absolute 
calibration and of the standard deviation of the fluxes in each 
radial bin. Typical errors are of the order of 10-30%. However, 
uncertainties for the dust mass are mainly due to the choice of 
the model parameters and can vary by a factor two. 

In the next Section, we report the total mass of H2 obtained 
with the different assumptions on Xqo (TabHJl. 

3.4. The Xco factor 

To convert '-CO observation to measure the amount of molecu- 
lar hydrogen we assume a CO-to-H2 conversion factor 



Xco =A^(H2)/ I(CO)dv 



(iLebrun et al.lll983h . In our Galaxy Xco- (1-8 ± 0.3) x 10^" 
cm~^ (K km s"')"' with excursions of up to a factor of ~ 2 
over this value, particularly at high latitudes ("Dame et"ar 1l200Th . 
LShettyetal. (2011) and Leroy et al. (2009) showed that Xco 
can vary even more, especially in low metallicity environments 
where it can reach Xco ~ lOOXca/. 



As derived by analysis of obs ervations by Bolatto et al. I 
(2008) and by theoretical models dGlover & Mac Lowl l20TTr 
this factor can vary due to effects of abundance, excitation, 
optical depth, and cloud structure averaged over a large area. 
In particular, the abundance of the heav y elements has an 
important impact on t he value of Xco dlsrael et al. I 119861: 
iMalonev & Black I 11988 ). Metallicity affects cloud structure 
both directly, as smaller abundances of C and O translate into 
lower abundance of CO, and indirectly, as a lower dust-to-gas 
mass ratio diminishes the H2 formation rates and the shielding 
of molecular gas from th e photo-dissociation b y ultraviolet radi- 
ation. As summarized bv lBolattoetXI (120081) . there are several 
observational calibrations of Xco with the metallicity, Z, in the 
literature, showing a range of behaviors; most of t he cahbrations 
find an increasing Xco with de creasing Z (e.g.. 'Wilson I I 19951 : 
Barone et al. 2000; Is rael et al.ll20 03: Strong et al. 2004)._^ 

The dep endence varies from Xco ~ Z' ^-^ (Israel 200"o') to 
Xco- Z ' (IWilsonlll995t lAri moto et al. |[T996: Boselli et afl 
l2002h . while^olatto et al. (200^) do not find any measurable 
trend in the range from 8 < 12 + log(0/H)< 8.8 on the scales 
of the individual CO-bright giant molecular c louds, but they did 
not argue for a Galactic Xco on large scales. IStrong & Mattox I 
(1996) also found a constant Xco with metallicity. However, few 
works were dedicated to the metallicity depend ence of Xco in 
super-solar regions, as for example the study of Arimot o et al. I 
(1996) who investigated Xco also in two metal-rich galaxies, 
namely M31 and M51. The spiral galaxies in Virgo are thus 
of particular importance to verify the assumed dependence on 
metallicity of Xco in metal-rich environments. To evaluate the 
conversion factor Xco at different metallicities and relate it to 
the dust-to-gas ratio, we compared three different assumptions: 

i) Uniform: we used the standard Galactic value from 
iBolatto etal.l (l2008h Xco- 1-8 x 10-" cm ^ (K km s ') ' 
(hereafter Galactic Xco ) throughout the radial range. 

ii) Linear metallicity dependence: w e used a linea r metal - 
licity dependence as found by lArimoto et al. I (Il996l) : 
iBoselU eTaTI (|2002|). This is similar to the -0.67 slope de- 
rived bv IWilsonl^995). W e adopted the analytic calibra- 
tion given by Wilson ( 1995), bu t rely on the m ost recent 
oxygen solar abundance of Asplund et al. I (l2009b instead of 
iGrevesse & Anders! (Il989l) . Assuming that at the new so- 
lar oxygen abundance, 12-H log(0/H)=8.69, Xco - 1.8 x lO^'' 
cm"^ (K km s"')"', the .Boselli et al. . (.200Z) 's relationship 
becomes 



log(Xco) = -1.0 X (12 -H log(0/H)) + 8.94 



(1) 



where Xco is expressed in unit of 10-" cm"^ (K km s"')"'. 
Hi) Super-linear metallicity dependence: we considered a super- 
linear dependence on metallicity as found by Israel (200(%. 
Re-scaling to the Asplu nd et al. I ([2009) solar oxygen abun- 
dance, the Israel . (.2000.) relationship becomes 

log(Xco) = -2.5 X (12 + log(0/H)) + 21.72 (2) 

where Xco is again expressed in unit of 10^". 

4. The method 

We started from the hypothesis that dust-to-gas mass ratio and 
O/H are tracing the same quantity, i.e. the metal abundance rel- 
ative to hydrogen. Thus we expected that the slope of the two 
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Fig. 4. Radial variation of position angle and inclination in NGC4254, NGC4321, and NGC4501 from moment-1 H I maps. The 
continuous lines are the average in the whole radial range, whereas the dashed line is the average in the region Rq < Rgpt for 
NGC4254. 



Table 1. Parameters of the sample galaxies. 



Name 


Inclination 


Position angle Do,,/' 




12+log(O/H)K^04 


Ref. 






(arcmin) 


(dex/Rzs) 


at Rg=0.7R25 


(7) 


(1) 


(2) 


(3) (4) 


(5) 


(6) 




NGC4254 


42° (53°)* [42°]" 


67° (61.6°)* [68°]" 5.6 


-0.42+0.06 


8.84+0.02 


MIO'^ 


NGC4303 


30°*^ 


135°'' 6.5 


-0.42±0.02* 


9.24±0.02 


S96" 


NGC4321 


33° [26°y 


154.4°[155°]/ 7.6 


-0.35+0.13 


9.04±0.04 


MIG"^ 


NGC4501 


63°[58°F 


141.4°[140°]'? 7.1 


-0.07±0.16* 


9.30±0.21 


S96" 



Notes, a D„,„ is 2XR25; b Average values for Rq <R25 and, between brackets, throughout the whole radial range, c M10='Moustakas et al. ' (2O10|) 
and S96=Skillman et al. ( 1996) ; d Phookun et al. ( 1993); e Koda & Sofue ( 2006); / Canzian & Allen ( 1997); g Kenney & Young 1 1988). $ 
We re computed the gradient slope of NGC4303 and NGC4501 using data from Skillman et al. ( 1996), converting them with lKewlev & Ellison I 
( I2OO81) formula to the KK04 scale, obtaining dilferent values with respect to those quoted in their paper. 



radial gradients agree. Within this hypothesis, we attempted to 
constrain the value of Xqo, and eventually its dependence on 
metallicity, using the information contained in the slope of the 
metallicity gradient. In fact, due to the large discrepancy on 
the zero-point of differe nt calibrations as discussed by, e.g., 
iKewlev & EUisonl (l2008h . it is difficult to place these galaxies 
on an absolute metallicity scale using only inferences from HII 
optical spectroscopy. 

Here we summarize the steps, while the details are given in 
the following sections. 

I We converted the dust-to-gas mass ratio in oxygen abun- 
dances with the relationship of Draine et al. (2007) (see 



Eq. O (i Draine et al. Il2007i see Eq. |3]l, thus obtaining an ap- 
proximate zero-point for the metallicity scale. 

II We computed the radial profile of the dust-to-gas mass ratio 
with the Galactic Xqo ■ 

III Once fixed the approximate abundance scale, we compared 
the radial profiles of the dust-to-gas mass ratio with the 
super-linear and linear metallicity dependences of Xqo ■ 



4.1. I. The dust-to-gas ratio and metallicity scales 

If we consider that the abundances of all heavy elements are pro- 
portional to the oxygen abundance and that all heavy elements 
condensed to form dust in the same way as in the MW, then the 
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Fig. 5. Top panels: radial profiles of dust emission at 100, 160, 250, 350, 500 ^lm (from top to bottom) for NGC4254, NGC4303, 
NGC4321, NGC4501. Bottom panels: radial profiles of dust mass(green solid line), HI (red long dashed line), H2 (blue dot-dashed 
line) surface density using a unifonn conversion factor Xco~ (1.8 + 0.3) x lO^'' cm"^, and the O/H radial gradient (dotted line, 
scaled of a factor 10.5 to match the logarithmic scale of the other quantities). 



dust-to-gas mass ratio scales proportionally to the oxygen abun- 
dance (Draine et al. 2007) CDraine et al. ..2007.) 



0.007 



(O/H) 
(0/H)o 



(3) 



where 0.007 is the dust-to-hydrogen ratio of the MW at solar 
radius, ~0.7R25 with optical radius of our Galaxy of 12 kpc 
dReshetnikov 2000). In eql3l we use the observed dust-to-gas ra- 
tio 0.0073 estimated from observed depletions in the solar neigh- 
borhood in stead of the value from dust models, 0.010, about 
40% larger (iDraine et al. I2007h. and the solar ox ygen abundance 
I2-Hlog(0/H')^^8.69 from lAsplund et al. I dlool . 




A similar relationship was obtained also b y iBoseUi et al. 
using the data available fo r the M W dSodroski et al. 



the LMC jKoomneefl [1981 . and the SMC 



(iBouchet et anil985h . 

We converted the dust-to-gas ratio, obtained with the 
Galactic Xqo at Rg=0.7 R25, into 12-H log(0/H) with Eq.[3] 

Before comparing the metallicity derived from dust-to-gas 
ratio with that derived from nebular oxygen abundances, we es- 
timated the effect of the assumed values of the mass emissiv- 
ity coefficient and of the Xco on the metallicity derived from 
Eq. [3] Because the determinations of the total gas mass take 
place in regions where H2 is not negligible, the dust-to-gas ratio 
depends indeed critically on what is assumed about Xqo ■ First, 



we checked the value of the dust-to-gas ratio we would obtain at 
Rc=0.7 R25 adopting different Xqo ■ using Xqo ~ 4 x 10^" cm"^ 
(K km s"')"' we would have a dust-to-gas ratio ~0.2 dex lower 
and using Xco ~ 0.5 x 10^" cm"^ (K km s"')"' we would obtain 
a dust-to-gas ratio ~0.2 dex higher. Thus an uncertainty of +0.2 
dex is associated to the choice of Xqo ■ 

Then, we verified how the emission coefficient could affect 
our discussion. In the most unfavorable conditions we could 
overestimate (or underestimate) the dust mass by a factor two, 
coiTesponding to a variation of O/H derived from Eqi3]of +0.3 
dex. 

In addition, we have that a typical error on the dust-to-gas 
mass ratio of ~30-50% (including the flux calibration) trans- 
lates into errors of ~0. 10-0. 15 dex in O/H. Considering that 
these sources of uncertainty are independent we combine them 
in quadrature, obtaining a total uncertainty of ~0.35 dex. 

We compared with 12+log(0/H) obtained at the same 
galactocentric distance with thre e cahbrations: KK04 , 
Pett ini & Pagel] (|2004 (PP04), and Pilvugin &TTiuanJ (l2005h 
(P05). The results for NGC4254, NGC4303, NGC4321, and 
NGC4501 are shown in Table |2] 

From the comparison of the metallicity obtained for dust- 
to-gas ratio with s everal metallicity calibratio ns we found that 
the calibrations of iPilvugin & Thuani ( |2005|) show the largest 
discrepancy with the metallicity derived from dust also if we 
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Table 2. Oxygen abundance from dust-to-gas mass ratio at 
O.7R25. 



Name 


^f^diist-lo-gas 




0/Hppo4 






0.7R25 








(1) 


(2) 


(3) 


(4) 


(5) 


NGC4254 


8.84 


8.95 


8.65 


8.35 


NGC4303 


9.24 


9.24 


8.94 


8.64 


NGC4321 


9.04 


9.10 


8.80 


8.40 


NGC4501 


9.04 


9.30 


8.90 


8.70 



consider all possible sources of errors, i.e., uncertainty on the 
emission coefficient, Xco, flux calibrations, and t end to give 
lower O/H values. The calibrations of KK04 and Pet tini & Pagel I 
(1200 4 are, on the other hand, both consistent within the errors 
(-0.35 dex 0/Hrf„ 

st-to-i;as snd ~0.1 dex for O/H from nebular 
abundances) with our findings from the dust-to-gas ratio. We 
conclude that, at the equivalent solar radius (0.7 R25), the four 
galaxies in analysis are more metal rich than the MW at the same 
radius (the solar oxygen abundance is 8.69), thus a perfect envi- 
ronment to test if any dependence of Xco with the metallicity is 
in place at high metallicity. 

Our resu lt is only in apparent disagreement with 
iDraine et al. I ( l2007h . Their relationship between dust-to- 
gas and metallicity is derived theoretically assuming that the 
interstellar abundances of all heavy elements were proportional 
to O/H, and their fraction scales as in the MW, and then it 
is compared with the oxygen abundances obtained with the 
calibration of Pilvugin & Thuan (2005) and the dust-to-gas 
mass ratio adopting Xco~ 4 x 10^" cm"- (K km s"')"'. They 
found a good agreement with the theoretical relationship within 
a factor of 2. In our analysis we found a better agreement with 
KK04 oxygen abundances since we are using a diff'erent Xco 
factor and a different value of dust-to-gas in the MW at solar 
radius: the P05 oxygen abundances are indeed on average 2-4 
times lower than the KK04's ones, but this is compensated by 
a Zco factor two times higher, and the assumed dust-to-gas in 
the MW for the dust models, ~ 1 .5 higher than the one we used, 
which comes from observations. 

With our choice of a Galactic Xco~ 1.8 x 10^° cm"^ 
(K km s"') a dust-to-gas ratio in the solar neighborhood 
of 0.0073, and the relation in eql3] we selected the metallic- 
ity scale of KK04, which best agrees with the scale derived 
from the dust-to-gas ratio, and convert the abundances of our 
galaxies to a common scale. For N GC4254 and NGC4321 
we used the abundance determined by iMoustakas et al. I (l2010l) 
(based on literature spe ctroscopy calibrated with the formula of 
iKobulnickv & Kewlev' ^04). For NGC4303 and NGC4501 we 
used the original spectroscopy to calculate oxygen abundance on 
the abundance scale of KK04. 

With the calibration of KK04, the abundances at O.7R25 are 
all super-solar ranging from ~8.95 to ~9.3, as shown in Table|2] 
The central oxygen abundances would range from 9.3 to 9.5, but 
we have to consider that they are an extrapolation of the metal- 
licity gradient up to the galaxy centre and not a real measure- 
ment. In addition, the galaxy disk does not extend up to the cen- 
tral regions where the bulge is present. These central abundances 
are comparable with recent results obtained from direct electron 
temperature measurement in the HII regions of M81, where the 



extrapolation of the metallicity gradie nt up to the galactic centr e 
gives an oxygen abundance of 9.37 (Stanghel lini et al.1l20IOl) . 
and in our Galaxy wh ere the central O/H from the gradient of 
iRudoIph et aTl (l2006l) is 9.2 (determination of O/H wifli optical 
spectroscopy). 

4.2. II. The dust-to-gas mass gradient with a standard Xco 

In Fig.|6]we show the comparison between the dust-to-gas mass 
ratio gradients, calculated with the Galactic conversion factor 
(red cui-ves), and the oxygen abundance gradients (magenta cir- 
cles and solid black lines) calibrated with KK04 and converted 
in the dust-to-gas mass scale with the relationship discussed in 
Sect. 10 

If we consider a constant Galactic Xco, the slopes of the 
dust-to-gas and metallicity gradients are only in marginal agree- 
ment. In particular, the dust-to-gas gradients are flatter, with a 
slight positive slope at large galactocentric radii Rc > O.7R25. 
ForNGC4254 and NGC4321 the comparison between the metal- 
licity and dust-to-gas mass ratio gradients is straightforward be- 
cause their metallicity gradients are well determined in the same 
radial regions where the dust-to-gas map is obtained (Fig. |6] left 
panels). For NGC4254 and NGC4321 the dust-to-gas mass ratio 
gradient are flatter than the metallicity radial distribution, and 
they are constant within the errors, with a slightly increasing be- 
haviour in the outer regions. The same is true also for NGC450 1 , 
while for NGC4303 dust-to-gas mass ratio gradient is positive. 
For these two galaxies we are planning new spectroscopic ob- 
servations of HII regions in the inner regions to have a better 
overlap between dust-to-gas and metallicity. 

4.3. III. Can we constrain Xco ? 

We used Eqs. 1 and 2 along with the abundance gradients given 
in Table [T] to calculate the molecular gas surface densities, and 
then the dust-to-gas mass ratios. 

The constant Galactic value of Xco (dot-dashed red cui-ves 
in Fig. |6l) produces in most cases flat dust-to-gas mass ratio 
as discussed in Sect. 14.21 The super-linear dependence (dashed 
blue curves) results in a dust-to-gas radial gradient much steeper 
than the oxygen gradient, especially in the inner regions. A lin- 
ear fit of the gradient in the logarithm of the dust-to-gas ratio 
[d log(dust/gas)/dRG] (Table O suggests that, within ~0.7 R25, 
the linear metallicity dependence is able to reproduce, within 
the errors, the same radial slope of the metallicity gradient for 
NGC4254 and NGC4321. A flattening of the dust-to-gas mass 
ratio gradient is appreciable in the outer regions, where the HI 
component is dominant. For NGC4303 a super-linear metallic- 
ity dependence of Xco is instead necessary to match within the 
errors the two gradients. For NGC4501 an Xco value between 
the Galactic one and the one derived from the Unear metalhcity 
dependence is required. 

From Fig|5]we can figure out what is happening: for a fixed 
Galactic Xco the dust-to-gas ratio would be constant, as in 
shown by dot-dashed red curves in Fig|6] Once a metallicity de- 
pendence of the type described in Eqs. 1 and 2 is introduced in 
Xco , that depresses the H2 in the centers (which have super-solar 
metallicity) and makes the total gas to look more like the atomic 
component, which generally has a much more constant distribu- 
tion with Rc than the dust. Hence, when adopting the metallicity 
dependences of Xco , we favor the outwardly decreasing gradi- 
ents in dust-to-gas ratio. 
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Fig. 6. Dust-to-gas mass ratio and metallicity radial profiles of NGC4254, NGC4303, NGC4321, and NGC4501 with the oxygen 
abundance calibrated with KK04. The conversion factors Xqo are: with the Galactic conversion factor (dot-dashed red curve), 
with the linear metallicity dependence as in Eq. 1 (dotted green curve), with the superlinear dependence on metallicity as in Eq. 2 
(dashed blue curve). The magenta circles are the oxygen abundance converted to the dust-to-gas scale with Eq.3. For NGC4254 
and NGC4321 we us ed the recalculation by Moustakas et al. (20 10) , for NGC4303 and NGC4501 the oxygen abundances by 
ISkillman et al. 1 (fT996l) converted to the KK04 scale. The continuous black lines show the fit to the abundance data. 



Table 3. Slopes of the O/H and dust-to-gas gradients within 
0.7R25. 



Name 


slopeo/H 


slopeo,,, 


Sl0pe;,„pa,- 


S 1 OpC , perl ill ea r 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC4254 


-0.42±0.06 


0.07±0.06 


-0.55+0.07 


-1.17+0.06 


NGC4303 


-0.42±0.02 


1.11±0.15 


0.27+0.13 


-0.44+0.11 


NGC4321 


-0.35±0.13 


0.09+0.09 


-0.52+0.08 


-1.13+0.09 


NGC4501 


-0.07±0.16 


0.07+0.07 


-0.45+0.07 


-1.18+0.12 



of the galaxies, while the outer regions, where the atomic gas 
components is dominant, are less affected by the choice of Xqo ■ 
From Fig. |7] where we show the radial dependence of Xco ob- 
tained with Eq. 1 (linear dependence) for NGC4254 (green cir- 
cles), NGC4321 (yellow triangles), NGC4303 (red diamonds), 
and NGC4501 (magenta squares), it can be seen that no larger 
radial variations of Xco are required. To reproduce a decreasing 
dust-to-gas mass ratio gradient in the radial regions where the 
molecular gas dominates the gaseous component, we need for 
the sample galaxies an Xco ~ 1/3 of the Galactic value, slightly 
increasing with radius up to ~2/3 the Galactic value. These re- 
sults agree with Xc o obtained in th e metal rich HII regions of 
M31 andMSl (e.g.. lArimoto etal.lll996h . 



Given that these galaxies are dominated by the molecular 
gas component in their central regions, one could think that the 
result we found is, in a certain sense, expected: we use a Xco 
that depends on 1/Z and we obtain a dust-to-gas mass ratio de- 
creasing as Z. However, the radial dependence of Xco has only 
a secondary impact on our results: to obtain a decreasing dust- 
to-gas mass ratio it is necessary only to have a conversion factor 
of the order of ~l/3-l/2 the Galactic value in the central regions 



To check this, we have used a constant Xco ■ 0.5 xlO^" cm ^ 
(K km s-i)-' for NGC4254, O.lxlO^" for NGC4303, O.QxlO^o 
for NGC4321, and 0.7x10^" for NGC4501. These values are 
chosen to produce the best agreement between O/H and dust- 
to-gas gradients for each galaxy of our sample. They are not re- 
lated to the Xco values derived from the metallicity dependence 
in Eq.l and shown in Fig|2l Note that the case of NGC4303 
seems to be quite extreme: the Xco necessary to reproduce with 
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Fig. 7. Radial variation of Xqo in the four galaxies with the 
metallicity computed with the calibration KK04 and the depen- 
dence in Eq. [1] NGC4254 (green circles), NGC4303 (red dia- 
monds), NGC4321 (yellow triangles), and NGC4501 (magenta 
squares). 



the dust-to-gas ratio the behaviour of the O/H gradient is much 
smaller than 1/3 Galactic Xqo , being 0. 1 x 10^°. 

We have compared these gradients also with the dust-to-gas 
mass ratio computed with constant "Galactic" Xco =1-8x10^" 



cm^^ (Bolatto et al. 2008) and Xco =4.0x10^" cm^^ (Draine et 
al. 2007). The results are shown in Fig|8] with the low Xco, the 
negative slope of the metallicity gradient is reproduced by the 
dust-to-gas mass ratio gradient. We note also that the outer re- 
gions are marginally affected by the choice of Xco , and they 
are still consistent with a "standard" Xco values for normal disk 
galaxies, while the inner regions need a lower value to repro- 
duce a decreasing dust-to-gas mass ratio. This is in agreement 
with the evidence f o und f or a lower Xcn in galaxy centers. 



have only a minor impact on the radial profile of the dust-to- 
gas mass ratio. Unlike NGC2403, in the galaxies of our sam- 
ple the molecular gas dominates the ISM within ~0.5 R25, and 
thus a va riation of Xco has a significant effect on the dust-to- 
gas ratio. iBendo et al. I ( 1201 Oa f) found with both the constant and 
metallicity dependent Xco ^ good agreement between the O/H 
gradient (-0.084±0.009 dex kpc"') and the dust-to-gas gradient 
(-0.093+0.002 dex kpc"'), with the dust-to-gas slope sHghtly 
steeper. They explain that variations in the ratio of oxygen to the 
other constituents of dust might originate the difference in the 
gradients. For example, carbon and oxygen are both produced 
by supemovae, but carbon is also generated by the nucleosyn- 
thesis of low and intermedia te mass stars; thus C /O is expected 
to vary radially, as shown bv lGarnett et al.l(ll999b . 

iThomas et al. Idlooi) obtained a similar result in their study 
of the radial distribution of dust and gas in several nearby galax- 
ies for which they had neutral hydrogen and 850 fim images 
available: their radial profiles of dust-to-gas mass ratio were 
flatter than the oxygen abundance gradients from HII regions. 
However in their galaxies CO observations were not available, 
so they were neglecting an important potentially dominant com- 
ponent in the inner regions. 

A similar result was obtained also bylM unoz-Mateos et al. I 
(l2009l) in the galaxies belonging to the SINGS sample. They 
have studied the radial variation of the dust-to-gas ratio in a 
large sample of galaxies of different morphological types. They 
found that the dust-to-gas mass ratio decreases by an order of 
magnitude from the center to the edge of the optical disk of 
each galaxy, similarly to what happens for the gas-phase oxy- 
gen abundance. Their relation between dust-to-gas mass ratio 
and oxygen abundance can be described with a linear scaling law 
for 12 + log(0/H)>8.9, as we found for NGC4254, NGC4321, 
and NGC4501. At lower metallicities, they found that the dust- 
to-gas ratios are systematically below this simple relation. They 
explain this fact to be due to the large amounts of gas that has 
not yet undergone star formation activity and that reside in the 
outer regions of spiral disks. 



Sodroski et al. 



T9941) . Ilsrael et all (l2006h . Ilsraell (l2009h . 



I Watanabeetal, (2001). On flie oflier hand flie dust-to-gas mass ^ ^ q^^^^ evidence for a low Xco 



ratio obtained with Xco =1.8x10^" cm"^ and Zco =4.0x10 
cm"^ tend to be flatter or positive, especially with the larger 
value of Xco ■ Thus, even if it is difficult to prove a dependence of 
Xco on the metallicity due to uncertainties on the O/H and dust 
mass scales, the comparison between the oxygen and dust-to- 
mass mass ratio radial gradients in the four galaxies under anal- 
ysis allows us to put a constraint on their Xco conversion fac- 
tors, favoring lower values of Xco in the radial regions Rc <0.6- 
O.7R25. 

In Table |4] we report the total masses of HI, H2, and dust. 
We computed the dust masses integrating the radial profile up 
to O.7R25. We remind that this limit is due to the requirement of 
a S/N> 10 at 500 jum, necessary to limit the uncertainties due 
to background subtraction and avoid the artefacts caused by the 
high-pass filtering in the PACS data reduction. The comparison 
with the results of Davies et al. (2011), obtained with an inte- 
grated analysis of the galaxies, indicates that only 30% of the 
dus t mass resides i n the o uter-most regions, with R>R25. 

iBendo et al. I (1201 Oah performed a similar analysis on 
NGC2403, a late-type galaxy with a relatively low molecular 
fraction. Because in NGC2403 the atomic gas is the largest com- 
ponent of the ISM, differing assumptions about Xco (they used 
a constant value and a metallicity dependence as in our Eq. 2) 



In agreement with our result, a value of Xco lower than ~2x 
10^" cm"^ (K km s"')"' was a lso expected in NG C4254 and 
NGC4321 from the analysis of lEales etaD (l2010h . They cali- 
brated the mass-opacity coefficient (Ky) using their dust-to-gas 
maps. Using a standard Xco value they found a mass-opacity 
coefficient significantly lower than the value (0.41 m^ kg"') ob- 
tained by James et al. (2002). To solve the discrepancy, they 
suggested that using only SPIRE bands they are over-estimating 
the temperature of the dust; if the temperature of the dust in both 
NGC4254 and NGC4321 were actually ^10 K rather than 20 K, 
this would be enough to increase the value of /<•,, to that expected 
from the James et al. study. However, this is not the case since 
with our fit of PACS and SPIRE bands we found dust tem pera- 
tures in the range 17< T <25 K found bv lEales et al. I (l2OI0i) (see 
Fig. 3). 

The other poss i bility they suggested to reconcile their results 
with iJames et al. I (l2002h was to have values of Xco -6 times 
lower than the initial value they assumed, ~2x 10^" cm"^ K km 
s"')"'. Considering that we are using a mass-opacity coefficient 
which is about half of that used by James et al., thus, to have 
agreement with our results, they would need Xco -3 times lower 
than their standard value. This is consistent with our result; val- 
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Fig. 8. Dust-to-gas mass ratio and metallicity radial profiles with the conversion factor Xco 0.1-0.9x10^° cm-^ (Kkms"')-' (see 
text for the Xqo adopted for each galaxy-orange solid curves), and with constant Xco =1.8x10^" cnT^ (red dot-dashed curves) and 
Xco =4.0x10^" cm"^ (green dotted curves). The magenta circles are the oxygen abundance converted to the dust-to-gas scale with 
Eq.3. The continuous black lines show the fit to the abundance data. 



Table 4. HI, H2, and dust masses 



Name 


M(HI) 


M(H2)™,f. 


M(H2)linear 


-M(H2 )siiperlinear 


M(dust) 




10' Mo 


10' Mo 


10' Mo 


10' Mo 


10^ Mo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC4254 


4.4° 


5.4 


1.8 


0.22 


7.0 


NGC4321 


4.8° 


4.4 


1.3 


0.13 


8.6 


NGC4303 


3.3* 


4.8 


0.8 


0.06 


5.1 


NGC4501 


2.7° 


3.8 


1.0 


0.05 


7.4 



Notes. a lChungetal.l ( l2009ah 

b integration of the radial profile shown in Fig. 3 of ISkillman et al. 



ues of Xco in the range ~0.5-lx 10^" cm"^ K km s"')"' can 
reproduce both the correct value of the mass-opacity coefficient 
and the radial profiles of dust-to-gas mass ratio. 

5. Summary and conclusions 

1. We analyzed FIR observations obtained with PACS and 
SPIRE, together with CO and HI maps from the litera- 
ture, of four spiral galaxies in the Virgo cluster (NGC4254, 



NGC4303, NGC4321, and NGC4501). We derived the total 
mass of dust from SED fitting of the FIR images, and used 
HI moment- 1 maps to derive the geometric parameters of the 
galaxies, which are disturbed by tidal interactions. Finally, 
we extracted the radial profiles of atomic gas, molecular gas 
and dust, and compare them with oxygen abundance radial 
gradients compiled from the hterature and placed on a com- 
mon abundance scale. 
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To avoid the large uncertainties on the zero-points of the oxy- 
gen abundance derived with the bright-line methods, we used 
the dust-to-gas ratio obtained with the Galactic value of Xco 
to fix a lower limit to the oxygen abundance. This allowed us 
to discard a set of metallicity calibrations and to reduce the 
uncertainty on the oxygen gradient zero-point. We converted 
the literature 0/H to the scale of KK04. 
We studied the dependence of Xco on metallicity by com- 
paring the radial gradient of metal abundance in gas phase 
(oxygen abundance in HII regions) and in solid phase (dust 
in emission from the FIR versus the total gas-atomic and 
molecular), assuming that these two quantities decrease at 
the s ame rate with radiu s. We considered a constant Xco 
(e.g.. iBolatto et al." l2008l) . Xcqo^Z ' ("linear" dependence) 
re.g.. iBoselh etal. I '2002V and XcqocZ'^-^ (' 'super-linear" 



dependence) (e.g., [Israel 2000). A linear fit of the gradient in 
the logarithm of the dust-to-gas ratio [d log(dust/gas)/dRc] 
shows that, within ~QJ R25, values of Xco in the range 



~0.5-0.9x lO'^" cm ^ (K km s ') ' were able to reproduce 
negative gradients, similar to the O/H ones for NGC4254, 
NGC4321, and NGC4501, while NGC4303 needs an ex- 
tremely low Xco value ~0.1x 10-" cm"^ (K km s"')"'. For 
NGC4254, these Xco values can be obtained with a linear 
metallicity dependence and oxygen abundance calibrated by 
KK04. For NGC4303 a superlinear dependence is instead 
necessary to match within the errors the two gradients, while 
for NGC4501 and NGC4321 a Xco intermediate between 
the Galactic one and the linear metallicity dependence would 
be necessary (always with KK04 O/H calibration). We sug- 
gest that a Xco lower than the standard Galactic one (e.g., 
IBolatto et al. 2008) is necessary in these galaxies to obtain 
a decreasing dust-to-gas mass ratio similar to the O/H gra- 
dient of HII regions. These low Xco values are favored in 
the radial regions Re <0.6-0.7R25 where the molecular gas 
dominates, while the outer regions, where the atomic gas is 
the main component, are less affected by the choice of Xco , 
and thus we cannot put constraints on its value. A large sam- 
ple of galaxies with available metallicity ad dust-to-gas mass 
ratio gradients is necessary to confirm and strengthen these 
results. 
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